RD-A146  155 


SINC  METHOD  OF  SOLUTION  OF  SINGULAR  INTEGRAL  EQUATIONS 
<U>  UTAH  UNI V  SALT  LAKE  CITV  DEPT  OF  MATHEMATICS 
F  STENGER  ET  AL  1984  ARO-19297  8-MA  DAAG29-83-K-0012 

F/G  12/1 


UNCLASSIFIED 


\>  .W_\>  "i  ’ll  --»  ?.AiK»i£l 


MICROCOPY  RESOLUTION  TEST  CHART 
utiowh.  naiiv  or  itmomos -m»- * 


aaHsmmgniizgisggm 


REPORT  DOCUMENTATION  PAGE 


,T%  JTT-74 1 * 


1.  REPORT  NUMBER 

ARO  19297. 8-MA  _ 


4.  Tl T lE  (ltd  SvOtltlt) 

SINC  Method  of  Solution  of  Singular  Integral 
Equations 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  RECIPIENT'S  CATALOG  NUMBER 


S.  TYRE  OF  REPORT  •  PERIOD  COVERED 


Technical  Report 


S.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR*  •; 


Frank  Stenger,  David  Elliott 


S.  CONTRACT  OR  GRANT  NUMBER*  A) 

DAAG  83  K  0012 


ERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

University  of  Utah 
Department  of  Mathematics 
Salt  Lake  City,  Utah  84112 


CONTROLLING  OFFICE  NAME  AND  ADORESS 

U.  S.  Army  Research  Office 
Post  Office  Box  12211 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  S  BORK  UNIT  NUMBERS 


12.  REPORT  DATE 

/9 


IS.  NUMBER  OF  PAGES 


MONITORING  AGFNCY  N 


CmntroltlnS  Of  tic*)  IS.  SECURITY  CLASS,  (ol  thlm  capon) 

Unclassified 

Is*.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


distribution  statement  (of  Ml*  Rap  on; 


Approved  for  public  release;  distribution  unlimited. 


IT.  DISTRIBUTION  STATEMENT  (a  I  Pm  a&arracl  an  farad  In  Bloc*  20.  II  0llmrmn  Na  RaporfJ 


is  supplementary  notes 

The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position,  unless  so  designated  by  other  authorized 
documents. 


IS  KEY  PORDS  (Continue  an  raaaraa  a! da  II  nacaaaary  and  (Banff*?  Of  Sloe*  mmOar) 


rr>  A' 


.  J 


20.  ABSTRACT  (Continue  an  raaaraa  alBa  II  naeaaaa ry  and  Manlfly  Of  Hock  i 


J  The  use  of 
solution  of  the 
arbitrary  index 
the  equation  is 
solution,  based 
Convergence  of 
shown  that  the 


the  Whittaker  cardinal  (or  sine)  function  for  the  approximate 
complete,  one-dimensional,  singular  integral  equation  with 
is  discussed.  Meekly  -we-  shall  be  'concerned  with  the  case  when 
taken  over  the  arc  (-1,1).  An  indirect  method  of  approximate 
on  the  equivalent  Fredholm  integral  equation,  is  described, 
the  approximate  solutions  is  discussed  in  some  detail  and  it  is 
error  decays  exponentially. 


Mo  /iirt.t- a»/> 


Site  METHOD  or  SOLUTION  OF 
HNCULAt  IVTECKAL  EQUATIONS 

Frank  Scantier* 

Department  of  Matheoatlcs,  University  of  Utah, 
Salt  Lake  City,  Utah  8*112,  USA 


David  Elliott 

Department  of  Nathamatlca,  University  of  Tasmania, 
Hobart,  Tasmania,  Australia. 


Abstract 

The  use  of  the  Whittaker  cardinal  (or  sine) 
function  for  the  spproslmate  solution  of  the  complete, 
one-dimensional,  singular  integral  eqti.it Ion  vl th 
arbitrary  tndes  Is  discussed.  Nostly  ue  shall  be 
concerned  with  the  case  when  the  equstion  Is  taken 
over  the  are  (-1,1).  An  indirect  method  of  approximate 
solution,  based  on  the  equivalent  Fredholm  Integral 
equation.  Is  described.  Convergence  of  the  approximate 
solutions  Is  discussed  In  some  detail  and  It  la  shown 
that  the  error  decays  exponentially. 

I .  Int roduct Ion 


In  this  paper  we  propose  to  Investigate  the  use 
of  the  Whittaker  cardinal  function  (or  sine  function) 
for  the  approximate  solution  of  the  singular  Integral 
equation 

a(x>w(x)  ♦  ^  j  “7™;  *  |  K(x.t)w(t)dt  •  f(x),  (1.1) 

for  x  i  f,  wti.tr  is  an  open  an  which  mostly.  In 
this  paper,  wt  take  to  br  (-1,11.  The  Turn  lions  a.  b, 
K  and  f  of  <1.11  are  assumed  to  be  given  on  V  and  It  la 
required  to  find  w,  or  approximations  to  w.  The  first 
integral  appearing  In  (l.l)  is  a  Cauchy  principal  value 
integral  defined  hv 


1  la  [  wft)dt. 
»  FT  *  "  * 


(1.2) 


whvre  r<  Is  that  part  of  T  cut  out  bv  a  circle  of 
radius  t  with  erncre  at  a,  provided  that  the  Halt 
exists. 


The  theory  of  equation  (l.ll  la  to  be  found,  for 
example.  In  the  book  by  Nuskhcllvhi  1 1 1  |8|  and, 
following  him,  we  shall  look  for  solutions  w  of  (1.1) 
which  are  absolutely  lntegrable  over  (—1,1).  Let  os 
define 


both  2  and  l/Z  are  absolutely  lntegrable  over  (-1,1), 
see  Elliott  |4|.  As  discussed  In  |3),  rather  than 
aolvr  (1.1)  for  w.  It  la  computationally  more 
convenient  to  solve  for  the  function  F  aay  where 

w  -  (Z/r)F,  (1.7) 

the  function  r  being  defined  by 

r(x)  •  (a2(x)  ♦  b*(a))^*  (1.8) 

and  aastaed  to  be  strictly  positive  for  all  x  In  T  and 
Its  end  points.  We  replace  (l.A)  by  the  following 
equation  for  F, 

(aZ/r)F  ♦  b«l(ZF/r)  +/t(ZF/r)  -  f.  (1.9) 

It  la  well  known  (see,  for  exaaqtle  (6)  and  (*))  that 
thta  equation  can  be  "regularised"  so  that  It  la 
equivalent  to 

F  -  b^((l/rZ)*L(ZF/r))  ♦  (a/rZ)3C(ZF/r) 

•  (af/rZ)  -  bdf(f/rZ)e  bF^,. 

where  la  an  arbitrary  polynomial  of  degree  x-1, 

it  being  understood  that  (  r  0  when  c  s  0.  Now 

( 1 .10)  Is  a  Fredholm  Integral  equation  of  the  aecond 
kind  and  wr  shall  use  It  aa  the  starting  point  for  all 
methods  of  this  paper.  Thus  we  shall  be  describing  an 
"Indirect”  method  for  the  approxUwte  solution  of  (1.1). 

A  review  of  the  use  of  sine  functions  for  the 
approximate  solution  of  various  functional  equations 
has  been  given  by  Stenger  19).  and  frequent  reference 
to  the  results  of  |9)  will  be  made  throughout  this 
paper.  However,  in  12,  we  derive  a  new  class  of  sine 
approximations  which  la  particularly  suited  to  the 
solution  of  (1.1).  For  the  Interval  (-1,1),  If  we 
let  JD  denote  the  eye-shaped  region  containing  (-1,1), 
are  (2.39),  then  wr  shall  assume  that  F  la  analytic  In 
£>  but  of  the  class  Lip  in  the  closure  otjb  j  we 
write  I  i  I  UE>)  (see  Mflnltlon  2.1*)  and  In  i2  we 

(S 

derive  Interpoletlen  and  quadrature  formulae  for  such 
functions. 


(*w)(i)  •  J  .  (*w)(x)  •  jn(,.t)v(«>dt  .  (1.3) 

r  r 

then  we  far.  rewrite  (1.1)  as 

sw  4  bjt  w  4  A  w  •  f .  (I.*) 

Suppose  that  (l.l)  has  Index  *  which  mav  be  positive, 
iw-grti  Ive  or  «cru.  (For  the  calculation  of  •  given  a 
and  b.  see  Dow  and  F.lllott  ())).  If  we  define 

C(x>  •  (a(x)  -  lb (x) ) / (a(x)  *  lb(x))  (1.5) 

then  the  fundamental  function  Z  Is  defined  by 

Z(x)  •  eap{-(l/2)(*log  fi)(a)J  .  (I.*) 


In  13  we  first  consider  some  of  the  properties 
of  the  Integral  operators  arising  in  equation  (t  .10) 
and  then  we  consider  Its  approximate  solution.  This 
ta  based  on  the  Calerkln  method  but,  as  we  shall  show 
In  I). 5,  the  sine  approximating  basis  that  we  use 
reduces  the  Calerkln  echcme  to  a  Ny Strom  achame  so 
that  It  Is  also  a  collocation  method.  The  details 
of  the  approximation  are  given  In  13.2  and  13.3. 

In  13.*  we  consider  the  convergence  of  the  approximate 
solution  to  tbs  exact  solution.  In  our  discretisation 
we  obtain  a  (ZH43)  term  linear  approximation  F„  aay  to 
F  which  satisfies 

eup  (Fix)  -T_(x)|  -  0(tfexpl-(ldoM>*&  (l.ll) 
-l«x<l  " 


where  we  choose  appropriate  branches  of  log  C  so  that 
it  Is  continuous  on  (-1,1)  and  furthermore  aurh  that 
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as  !•».  It  is,  of  course,  this  exponential  (rather 
than  algebraic)  decay  of  the  error  with  H  which  makes 
the  uee  of  elnc  function  methods  in  maerical  analysis 
eo  attractive.  Although  we  shall  not  do  so  here.  It 
nay  be  shown  that  the  rate  of  decay  on  the  right  band 
aide  of  (l.ll)  is  optimal  (ooe  (2))  in  the  sense  tbit 


» 


there  1*  no  baal*  lu.M),  k  *  (-N-I)  (l)(MI), 

I  ■  1,2,3,...  such  iMt 
Ml 

•up  | P(m)  -  T  e.(F)uLM(»)| 

-1<«<1  k— *1-1 

-  0(P<N.N‘I>'*I,<-3*|S>).  (IU2) 

M  *■••,  whcrt  p  If  a  fixed  polynomial  in  N  and  M**, 
mod  where  y  >  (adn)  \ 

Final  1 y ,  In  $4.  we  consider  •  particular  cioaplc. 
This  paper  la  latendvd  to  provide  the  fimmJ.it ion  for  the 
application  of  alnc  act  hod*  to  ilte  approxlaatc  solution 
of  singular  Integral  equations.  Much  rrsilm  to  be 
done  and.  In  particular,  one  Might  scntlon  tlie  study  of 
direct  artltoJn  for  approslaatr  solutions,  eatenalonR  to 
other  arcs  (only  briefly  Mentioned  here)  and  closed 
contours,  and  finally  the  application  of  these  Methods 
to  syateMS  of  singular  Integral  equations. 

2 .  Integration  Fonsulaa  for  Evalu.it 
Cauchy  Principal  Value  Integrals. 

The  definitions,  notations  and  results  of  thla 
section  are  taper tant  to  the  rest  of  the  paper.  We 
•hall  derive  two  faalllea  of  foraulae  for  a  very 
general  contour  T.  These  two  faalllea  reduce  to  a 
single  faally  In  the  case  when  T  ■  I.  The  laportant 
case  of  I  ■  (-1,1)  ta  then  given  special  consideration. 

2.1.  The  Doaalnf)^  and  Approalaat Ing  Functions  on  R. 

Definition  2.1.  let  R  denote  the  real  line, 

•  •  <  —  ,«*)  let  C  •  {x  •  *  ♦  ly:  k  •  t,  y  i  R)  and  let 
Z  •  {k:  k  •  0,  SI,  i2,  ...}.  Let  d  and  h  denote 
positive  n  iBbria  and  let  ua  define 


-  P(»k) 

'  *  J.  F55? 


F(«)e*pU(*4(t)/h)»qnlM 

ein|»4(t)/h) 


(a)) (da  . 


|r  <  f  :  | far  1  <  dj , 

(2.1) 

•ln].(a-kh)/h[ 

(M.-khl/hi’  •  *  *  z  * 

(2.2) 

1  -  eo»l*(H..*)/h)  .  , 

“F(kh-T»7hr  •  k  '  *' 

(2.3) 

2,2  The  More  General  Don. In  aft  and  the  Conformal  Map  t. 

Definition  2.7.  Lr(£  br  a  almply  connected  domain  In 
thr  complex  plane  C,  and  drnote  by  >£  the  boundary  of 
Jt)  .  Lot  a,b  (b  *  a)  ba  boundary  point*  of  eft  and  lot 
♦  be  a  conformal  map  of  fu  onto  (are  (2.1))  auch 
that  *(•)  -  -e,  *(b)  »  ».  Let  I  •  (*■  denote  the 
Inverse  map  and  act 

r  -  {4r(x>  :■<*},  (2.*) 

rk  •  *(Vh)  ,  b  r  *  .  (2.3) 

.Let  lit)  denote  the  family  of  all  function,  F  that  are 
analytic  In  X)  and  such  that 

■(F.fc)  •  [|r(f)di|  '  inf  f|F(t)dt|  (2.b) 

U>  C.*.C'*jc 

Theorem  2.3.  Let  S(h,h)  and  T(b,h)  he  defined  by  (2.2) 
and  (2.3)  re.pec lively  and  lat  *  >  I ( d.) .  Then,  for 

oil  x  c  r, 

;'<l  *  J,  rrv 


|  r(x)s(k,h).#(«)dx  -  hr(tk)/*,(xk) 

*  ,h^'-)k  I  r(«)e«pU(*»(r)/h)om(l,»(»))ldx. 

JL  ♦(*)  •  kh 

**  (2.*) 

and 

i  f  Tilrkr*)  -  jz  ^,T(k.h).*(x)  -  £ . 

r  h 

if (i )lcox|nf (x)/hl  -  exp|l(r»(t)/h)*xn(lm»(t))l dx . 
l*(x)-*(x) )ainlx*(i)/h) 

(2.10) 

Moreover,  If  the  left  hand  oldca  of  oquatlona  (2.7), 
(2.S),  (2.3)  and  (2.10)  are  denoted  by  •>,(>),  hj.  Oj, 
nt(i)  reapec lively  then 

|n,(x)|  i  M(r.«>)/(2»delnh(ad/h)).  a  t  T  ,  \ 

|n.|  <  exp(-xd/h)N(F,®  )/(2  etnb  (ad/h))  .  / 

1  V2.ll) 

Ih,!  -  h  exp(-ed/h)H(F.»)/(2*d)  ,  f 

|«i4l*)|  «  (l«exp(-iid/h))N(F.<©)/(2ed»lnh(»d/h))}) 
a  .  r. 

Frool .  See  Stenger  I») . 

Assumption  2.* ■  In  addition  to  F  r  S(b).  let  ua 
aaaumc  that  for  all  x  t  T, 

|F(x)/*'(x)  |  «  C  e*p(-o|b(*)|),  (2.12) 

idirre  r,a  ore  positive  ovmbers. 

Theorem  2.3.  Let  d(  H  denote  the  left  bond  side*  of 

aquations  (2. Ml),  I  •  1,2  and  *.  for  thr  esae  when  the 
Infinite  Bums  lk(I  are  replaced  by  finite  sums 

lot  T  aatlsfy  (2.12).  If  b  Is  ooloctad  by 

the  formula 

h  •  (xd/ah)1*  (2.13) 

than  tHrrt  exist  conataota  C.  which  are  Independent  of 
li  our h  that 

1*1,1  ’  ® ,,^eapj- (xdoN)^} ,  l  .  1,2,4.  (2.14) 

If  h  la  saloctod  by  the  formula 

h  •  (2*d/oll)^  (2.15) 

than  there  exlata  a  constant  CJt  auch  that 

l«2>J  *  Cjaapj-IZadaM)1*) .  (2.16) 

If  h  la  aelaetad  by  the  formula 

h  •  T/M1*  (2.17) 

Where  >  la  a  pealtlve  eoaatant,  than  there  ealat 
posit Ive  constant*  C*  and  <  such  that 

|«,  „|  *  C'Mp(-6nS,  1  •  l.J.3.4.  (2.U) 

Froof ■  See  (temper  | 9] . 


V] 


JD  i 3 1  r 1 but 1 on/  J 

Availability  Codes 
jAvali  and/or 
)ist  Special 


Mt  observe  chat  the  bound*  on  |n,(s)|  and  |n^(a)|  In 
Theorem  2.5  apply  uniformly  for  all  a  »  I .  moreover 

tha  Integral  (I/*)  j  (F( t )/(*(t >-9(x)))dt  a*  well  aa 

It*  approximation  via  (2.10)  approach**  trro  a*  « 
approach**  an  end-point  of  I".  Thl*  1*  not 
a*cra*arlly  true  for  the  new  approxlnat Ion  which  we 
•hall  derive  next. 

Let  ua  define  the  function* 

V(X)  5  aU.k.h.x)  .  .  <2.,,) 

t.  (*)  1  t(».k.h.<)  .  •  (2.20) 

*  ‘  'k 

where  we  intend  to  u*e  the  abbreviated  notation  when 
the  role*  of  1  and  h  are  underatood.  Froo  equation* 
(2.2),  (2.1).  (2.19)  and  (2.20)  It  follow*  that 

S(k.h).»(«  )-{“•}{  J  *  *•  ] 

M.  11  *  t.  1(2.21) 

T(l  hl.»f,  1  ./(l-(-l)l“,)/(*(k-t)).  If  k»l,J 

T(k.h).»(«t)  -|0  „  k  .  ,.  J 

and 

.  ,  (0.  If  k  »  I. 

V V  *  1  *'(*,).  If  k  •  I. 

.  ,  (2.22) 

,  .  ,  fh(!-(-l>  >/<*(»,-».)).  If  k  »  f. 

VV  1  0.  If  k  .  f.  k  1 

Theorem  2,».  let  F  >  ■(£)).  Then  for  .ill  x  /  f, 

x  »  a.b 

F(V 

t  (x)  ?  Fix)  -  l  rt77-.*k(x) 
k.  7  ’  1  k' 

.  *ln(x«(x)/h)  f  FO)dr _  , 

Jrl  J  (r-ilalnhUriyh]  * 

Cj  ‘  |  F<xtS(i(x)dx  -  hF(xk) 

r  . 

•  >h<-l )  f  F(x)expll(e«(t)/h)*gnlml(x)l 

2"  J*  *  "  *k  (2.2*) 

«><■•  •  i  j  W-  -  J,  ?&•.<•> 

.  _L  f  F(x)|co»|et(x)/hl-exp|l(vt(x)/)i)*gnl»d(«)]ld». 
2rl  j  (*-*)  aln  |*9(t)/h] 


It  ,(*) !  Nir.JD  ,x)/(2x*lnhfrd/h)),  (2.26) 

[•2!  -  C>/2-)e«p(-xd/h)h(F.®,rk),  (2.27) 

|f,(x)|  N(F,JD,*)|leexp(-»«l/l.))/<>«lnb(rJ/li)), 

(2.26) 

where  N(F,iD,x)  •  I  |F(z) /(r-x)dr | .  (2.29) 

W 

Proof .  To  prov*  (2.23)  loi  u»*  In  (2*7),  rrplirr  n 
by  I  and  F(«)  by  ((*< O -+(«))/<(•»)  |f(t)  to  give 


iiUltliMisii . 

9'(t)(t-x) 

*lnl*9(*)/hl 

2x1 


F(z)(6(z)-6(x))dz 
16(*)-#(t)  l(*-a)ala  {*♦(*) /i»l  ’ 


thl*  equation  holding  for  any  F  <  1(0)  and  all 
interior  point*  x  and  t  of  T.  Sotting  x  »  t  In  (2.10) 
give*  (2.21). 

Equation  (2.2*)  follow*  directly  from  (2.8)  If,  tn 
(2.6),  we  replace  F  by  *lF  and  uae  definition  (2.19) 
for  *k>  noting  that  iln|(t/h)(f(x)dih)]  • 

•  (-l)kaln(«*(x)/hl. 

To  prove  (2.25)  we  multiply  both  aide*  of  (2.30) 
by  t'(t)/(6(t)-4(a)).  Integrate  over  T,  and  uee  the 
ldentltle* 


(  £Lni!lZ!)igt  •  xexp{ (l«z/h)agn  Imz}, 

)g  '•* 

i  t  C,  ■  f  I  and 

f  «l«(nt/hjdt  .  ,ro.(lw/h)> 

'6  * 


It  •  »eo»(*x/h). 


x  «  6,  from  which  the  reault  follow*. 

Finally,  the  error  bound*  (2.26),  (2.27)  and 
(2.26)  follow  directly  by  bounding  the  contour 
integral*  in  (2.23),  (2.2*)  and  (2.25)  respectively .  f 

He  remark  that  under  our  aeaunptlons  on  F  the 
right  hand  *lde*  of  (2.23)  and  (2.25)  nay  became 
unbounded  e*  x  approaches  an  end  point  of  T.  Bence 
(2.23)  to  (2.25)  nuet  be  Interpreted  to  be  accurate  in 
the  *en*e  of  a  relative  error.  An  abaolute  bound  1* 
possible  l.e.  (2.23)  and  (2.25)  hold  for  all  x  on  T  If 
H(f.<0)  :  sup  H(F,£),x)  <  »,  nee  (2.29).  In  general 
x  t  F 

we  must  exercise  caution  In  the  evaluation  af  the 
approximating  sums  In  (2.23)  and  (2.25).  For  example. 
If  x  •  i|,  then  Cj(xt)  •  0  which  poee*  no  problem. 

On  the  other  hand  If  we  define 


(*F)(x) 


I  I'  Mt)dt 

w  J  t-*  * 


then  C*  F)(*f)  It  not  MctMirlly  itro  and  %rc  hive,  by 

<*•»>  k_, 

-hr  1 

(*F)(»t)S£  J  >  - — .  (2.3*) 

1  *  k*t  ♦  <*k)(V*t) 

The  evaluation  of  thl*  sun  to  within  a  relative  error 
of  (  c»n  be  carried  out  by  means  of  Algorithm  2.7  (see 
below).  The  approximation  of  F(x),  *  *  *,  by  (2.23) 
can  be  carried  out  by  a  similar  algorithm. 

Algorithm  2.7.  Evaluation  of  the  eua  In  (2.3*)  to 
within  ?  relative  error, 
k  •  I  -  1 

S,  •  2F(*k)/|9,Uk)(*k-*t)) 

«'  -  V  -  I;  H  •  |S.  | 

k  •  k  -  2< - 

T  -  2f ( *.  ) / 19’  (*.)(*.-*,) ) 

I.  •  I.  *  T 
0  •  V,  V  •  H,  H  >  |T| 

((’  ♦  V  ♦  H)/|$,  |  !  6 


I 


,r 


(<)  k  -  t  ♦  I 

Sj  -  2r(*h)/|*,(ikHik-*l)j 

0  -  1.  W  -  I.  V  -  |s2l 

k  •  k  ♦  2  q— - 

T  -  2F(rk)/|f(ik)(jk-it)l 

S2  -  S2  ♦  T 

H  -  V.  V  -  W.  W  -  |T| 

(U  ♦  V  ♦  W)/ |Sj|  :  6 

(«)  S  -  (h/w)(S,+S2) 

In  order  to  replace  the  infinite  sum*  In  Theorem 
2.6  by  finite  Bums,  we  consider  three  special  cases  of 
the  transformation  ♦  . 

Ea  .2 ,B .  r  ■  (-1,1),  In  this  case  choose 

♦“>  *  •°«(rH)  •  *k  *  ^  • k ' 2-  ? 

*  >(2.3i) 

&>  -  (i  f  C  :  |arg|j^jj|  <  d).  0  <  d  <  x.J 

E*.2.9.  r  •  (0,“).  Choose 

♦  (*>  •  lo|  «  ,  x  -  *kh  ,  k  •  I  .  i 

„  ,  f  <2. 3b) 

SU  •  |t  >  C  :  |arg  t |  <  d)  0  <  d  •  n .  J 

Although  thr  m-ipptng  (2.36)  suffices  for  anny  problem* 
over  ,  it  is  unsuitable  for  I  Itos**  lor  whl*  l»  llw* 

corresponding  functions  to  be  approx la.it til  are 
unbounded  in  the  reg  I (>ii £>  of  (2.3*').  Tin*  following 
transformation  is  suitable  In  the  case  whin  these 
functions  are  analytic  and  bounded  in  a  strip 
containing  the  interval  (0,**) . 

E* . 2 . 10.  T  •  (0,-).  Here  we  clioosc 

«(»)  *  logfalnht)  ,i  .  log|ekh.(l.e2k',)',l,k<z) 

„  .  1(2-37) 

^  •  1*  :  |arg(sinh  i)\  <  d),  0  <  d  <  n .  J 

Graphical  illustrations  of  the  regions  JD  of  (2.33)  to 

(2.37)  are  given  in  19). 

L— a  2.11.  Under  any  of  the  transformat  ions  (2.35)  to 

(2.37)  we  luive 

•up  |.  <*>/»'(*  >1  s  rh  ,  (2 . 38) 

a  t  r 

•up  |t  (*)/♦*(») |  <  »h  .  (2.39) 

•  i  r 

Proof.  To  prove  (2.1»)  we  ronxl.hr  ffrM  x  to  ho  In 
the  Intorval  l,k.|»*k,| I  •,,<l  then  *  »«•  ho  In  (ho 

remainder  of  T .  For  *  j  <  a  5  tk+),  wo  have 

|aln{ (n/h) |#(x)*kh) | |  «  («/h)|»(x)-kl>|  -  (2,40) 

Hence,  nn  thl*  Interval,  for  any  k  <  Z 


On  thr  lntrrval  l"\ l*k.i>*g.il*  **  ***  th*  Inequality 
|aln((v/h) Ir(x)-kh) ) I  <  1,  than  for  nay  k  a  l 

!£!  L  „„  h/» 


C,r|Tk-l)h.(k-rl)h]  l*'(*k>l,’,(Cl,l  * 
tec  Definition  (2.2). 

•  ly0  I .  iiiim  |  ,  i>  a  *2) 

Iron  (2.41),  for  aaeh  of  the  traneforantlonn  (2. IS)  to 

(2.37). 

Cenelder  now  (2.39).  On  the  interval 
l*k-l**k.l'  *"  h,v* 

|lk^  j  |h|l-coa|(x/h)(»(x)-kh)Hl 
I*  (*k) I  I  k*'(tk)(x-xk)  I 

if  ♦,(t)aln{(x/h)  (♦(t)-kh]  |dt  | 


if  ♦,(t)aln{(x/h)  |*(t) 
I  ♦TTxk) (»-*k) 


f'l*. 


,  for  k  r  Z. 


•k(x) 

♦  <x)-kl. 

♦  *(Tk) 

fWxk)(X-fch) 

Celt  x  'V 

w,*k-rV»)  k 

for  aaeh  of  the  tranaforaatlena  (2.3S)  to  (2.37). 


-'*k-r‘k' 

Nenre  as  in  (2.41)  we  have 

|tk(x)/*,(xk)|  i  a\  (2.43) 

Finally  on  rU*k.|**k^]J  »*  have 

|l  -  roa{(x/h) |q(x)-kh)} |  *  2  ao  that,  for  k  r  Z. 

|hl!l  I  <  (2h/w)  2  h  2 

U'(*k)|  !♦’ (ak» (i-*k)  |  *«*  •  (2-‘4) 

where  the  laat  Inequality  followa  by  proceeding  ••  in 
the  derivation  of  (2.42).  ( 

Aaeuaptlon  2.12.  Let  oa  ainae,  In  addition  to 
PcTfl D),  that  for  none  poaltlv.  conatanta  a  and  C 

|F(a)|  *  Ca^hi*)!  t  ,  ,  $  .  (2.45) 

Me  no to  that  (2.4S)  la  equivalent  to 

|P(x)  I  *  C  |1-*2|“,  a  .  (2.4b) 

for  the  ca*e  of  (2.3S);  to 

|P(r)|  *  C  |t/(H-t2)|“  ,  x  t&.  (2.47) 

for  the  race  (2.36);  and  to 

|P(*)|  *  C  |ta-I|°  ,  t,B.  (2.46) 

for  the  case  of  (2.37). 

Theorem  2.13.  If  Aaetaption  2.12  fa  aatlafiad  and  If 
6  la  any  one  of  the  tranaforaatlena  defined  In  (2. IS) 
to  (2.37),  then  by  taking  h  -  lad/aH)’,  there  eaiat 
conatanta  C,  and  CJ(  lodapandant  of  H,  Ouch  that  for 
alt  a  r  T 

N 

|F(x)  -  l  (F(tk)/6'(ak>)ak(*)|  t  C,HSn»p(-(i»dalt)’), 

k*'N  (2.49) 


|F(*)-(LF)(*> 


and 


N  K'lI 


k— a 


T  ,  k-  t 


(«) 


(2.50) 


Proof.  If  Assumption  (2.12)  la  satisfied  chrn  It 
follow*  let  each  of  the  transformation*  t  of  (2. IS)  to 
(2.37)  that  M(F.JP)  t  sup  -  mr,A>,a)  <  -,  aee  (2.2*). 
Then  for  (.(>),  J  •  1,3  defined  by  (2.26)  and  (2.2» 

It  follow*1  that 

|r j <■ )  |  i  Aeap  (-*d/h).  (2. SI) 

where  A  la  Independent  of  h.  Hence,  for  (2.69),  we 
havi*  on  uh  1  ng  (2.10)  that 
M 

|r(i)  -  )  rd  )(•.(!>/♦*(*  >)| 

k—N 


A  *-»d/h  +  ,h  c  ,r(  , 

|M*»  k 


Ar-’d/h  4  K.h  I  e-“kh  .  fro.  (2.66) 
k-N+1 


-  Ar‘,d/h  ♦  2Ceh.^(H+,,h/(l-*^h) 


|F(a)-(LF)(a>  -  f  |F(l.)-(LF)(»k)]S(k.h)4*(«)| 
kM 

*  C  N^aapf-fadodl)1*) ,  (2.52) 

|F(*)-(LF>(»)  -  \  |F(ik)-aF)(«k))(ak(a)/*,(ak))| 

k»-a 

<  CH>>*ap(~(*dUI)S.  (2.5d) 

lS{‘££41 . - 

N  .  , 

J  |F(»k)-(I.F)(r.  )JT(k, !.)•♦(*)  |c  (Veapf-fadoN)’), 
k-K  (2.5*) 


and 


I  F(t )dt 


t-a 


-*(LF>(a)  - 


k— N 


i, 

JIF(rk)-(LF)(*k)](tk(*)/*,(ik))|*cA*p(-(ado»)S. 


(2.60) 


Froof  ■  Thl*  la  a  consequence  of 
Theoren*  2.5  and  2.13.  I 


2.15  and 


<  A.-,,d/h  ♦  TCe^^/loh). 

Inequality  (2  .51)  now  follow*  on  choosing 
h  •  (nd/(uN))*.  The  proof  for  |t.(a)i  follow* 
olmllsrly,  using  (2.39).  I 

2.3  The  Special  Case  of  f  «  (-1.1).  In  this  cane  the 
transforaation  p.  the  region  &  and  the  point*  i.. 
k  <  2  are  defined  a*  In  (2.35).  Thr  region tf  fa 
bounded  by  two  circular  area  Intersecting  at  the  end 
points  !l  at  an  angle  d  aay  (aee  Figure  6.2  of  (9) > . 

If  F  >  Cl-I.ll  we  define  LF  by 
(LF) (f )  -  ( l-t)F(-l) /2  4  (l4f)F(l)/2.  (2.52) 

Then,  see  (2.33), 

(*I.F)(a)  -  (1/«)|F(-1)4F(1)1  -  (l/n)9(a)(LF)(a). 

(2.53) 

Definition  2.16_.  Let  «£*  be  drrined  a*  In  (2.35),  let 
0  e  j  v  |  _  mid  let  fi  (<t>)  denote  the  fnally  of  all 
fun'-floii*  tint  are  analytic  In  rft1  and  o(  cla**  Lip  In 
,  I  lie  i  I*. .urn  of  A’. 


kemnrk  2.17.  If  h  la  aalactod  by 

b  •  y/H*  (2.61) 

where  y  la  a  positive  constant,  than  thare  aalat 
positive  conatanta  C.  and  (,  Independent  of  M,  ouch 
that  thr  right  hand  tides  of  (2.57)  to  (2.6ft  nay  ha 
replaced  by  Claspin’1) . 

5.  The  Integral  Equation 
3.1  The  Operator  T. 

Following  the  dlocusslon  of  II  (oae  (1.10)), 
let  T  be  an  operator  defined  by 

Tf  •  •  b#(^*(T))}*  o  l) 

Let  u*  define  k(a,t)  by 

k(n,t>  •  K(*,t)Z(t)/r(t).  (3.2) 

AssuopMcn  3.1.  Let  £)  and  B _(•&)  be  defined  as  In 
IMinl'iTui  7.16  and  let  u*  assume  that: 


Leona  2.13.  Let  F  <  B  (A’),  then  the  function  C 
where  u 

C  •  F  -  LF  (2.56) 

Is  analytic  inland  satisfies  (2.66)  for  all  a  tJS . 

Proof .  We  hive 

|F(r>  -  F(l)  |  »  K|l-*|°  .  t  r  tF  . 

|F(»)  -  F(-l)|  «  k|l4* |“  .  t  ,  <£  . 

where  k  I*  the  Llpachltt  constant  of  F.  Now  if  it ti 
and  |l-x|  '  I,  then  |let|  ?  I  so  that  for  all  auch  a, 
«|l-t|a  s  It)  I-*2!'* .  Slol  lsrly  If  1 14«  |  '  1,  t*ien 
) I -a  1  t  1  and  the  1  ew  follow*.  f 

Theorem  2.16.  Let  *,  *k  and  A  be  defined  aa  In  (2.35L 
Let  F  <  Ba(E>),  let  h  be  chosen  no  th.it 

h  •  (*d/(*0>\  (2.5h) 

and  let#f(LF)  be  defined  by  (2.53).  Then  flu-re 
taints  a  constant  C,  Independent  of  K,  *u.  h  tint  for 
all  a  r  (-1,11, 


(2.55) 


(a)  frr  tone  flaed  a,  r  (-1,1), 

rl 

«  S  Jjk(*,,t>|dt  <  -  ;  (3.3) 

(b)  k(',t)  la  analytic  la  ©  for  sack  flaed 

I  (  (-1,1); 


(c)  k(',t)  satisfies  a  Llp#  condition,  l.e. 

|k(*,t)  -  k(y,t) |  *  C(t) |a-yl® 


(3.6) 


for  all  t  <  (-1,1),  for  all  a,y  c  fb  and  where 
C(t)  la  such  that 


T  *  C(t)dt  <  •  . 

3.2.  Let  Q  be  defined  by 
(OP) (a) 


t)P(*)dt. 


(3.5) 


(3.6) 


Then  Q  la  a  compart  operator  mapping  L  (-1.1)  into 


Proof .  tot  T  <  L  |-l,l|  bo  (Ivon,  and  let  us  set 

U F IL  •  sup  I F<«> | .  (1.7) 

Than  tlir  Inequality 

|k(>,  t)  |  *  |h(s| >t>  |  ♦  |h(a,t)  •  b(*j,t)  |  O.*) 

combined  with  (3.)),  (3.*)  and  (3.3)  yield* 

|(qr»(«)l  <  <6+*I°3)||pC  0  *> 

where 

i  I  f  2(  if  0  «  d  '•  n/2,  < «  1  n\ 

"  ,‘*y  *l2/*tnd- ,f  1,12  •  d  * «. 

y<  **>’ 

from  which  we  have  that 

| q ||  ■  6  ♦  nV  0-**» 

Moreover,  we  clearly  have  tha t _QF  is  analytic  in  £> 
and  furthermore,  for  all  x.y  r|)ve  have,  by  (3.4)  and 
(3.3)  that 


| (QF) (x)-(QF) (y) |  < 


<  ['  |k(x, t)-k(y,t ) |  |P(t)|dt 

*  Y|x-y|“|Fj,.  <3 


That  is,  QF  anp»  a  bounded  art  { F  c  L  (-1,1):  lrL « *1 
Into  a  family  of  functions  that  are  analytic  and 
uni  fore 1 y  bounded  in  and  of  class  Lip  lr»8>.  This 
proves  Leams  1.2.  #  a 

Let  3.3.  Let  a  and  b  be  in  B^(J&)  and  suppose  that 

|r|  •  !<«2*0S|  >  «  -  0  In  Si  .  (3.13) 

Then  the  operator  T,  see  (3.1), is  a  coap.nl  operator 
tapping  L* | - 1 » I  ]  into  Ba($  ) . 

Proof  ♦  Setting  f  •  QF,  it  suffices  In  view  of  Lrt 
1.2  to  show  that  the  operator  II  defined  by 

is  a  coapact  operator  tpping  B  (Jfc*)  Into 
Now,  see  (4),  the  solution  to  the  equation 

aZg/r  ♦  bjfdg/r)  •  f  <3. IS) 

having  Index  r  is  given  by 

*  •  *f  ♦  bp,-i  <3.U) 

where  is  an  arbitrary  polynowl.il  of  degree  r-1 

(P^l  I  0  if  •  s  0).  Also  under  the  given  conditions 
g  «  Llpft  1-1. 1 1  whenever  f  •  L i 1-1,1 1,  see  (5). 

In  order  to  show  that  Rf  *  B  ) ,  ler  P  be  an 
arbitrary  constant  In  tin*  range  -d  •  P  *  d,  and  let 
A(B)  denote  the  clrrular  arc 

*(?>  •  |t  <£  :  org(l+’.)/(l-C)  •  P).  (3.17) 

Let  u«  fix  t  »  JD  and  (  xurh  that  P  «  ar,(  Ur) /( 1-x) 
and  ronaldrr  (hr  function 

w(t>  •  ^  j  !  k  •  i  .  (3. Id) 


On  letrlng  «  ♦  C  *  A(6)  we  gat 

w(t)  •  b(t)|lh(C>  ♦  £  | 


Taking  P  •  -d  In  (3.IR)  wt  eee  that  w(t)  la  analytic  In 
O  .  Next,  with  (  an  arbitrary  point  of  ,  w*  oaa 
fro*  (3.19)  that  alaca  both  *(()  sod  b(L)h(()  arc 
analytic  In  Jb  ao  la 

wo  i  ^  j  CJ-iO) 

A  (6) 

Ve  nest  show  that  Rf  r  Lip  (Jb  ) .  To  this  and.  It 
oufflres  ro  show  that  the  Solution  g  co 


s(x)Z(»)g(») 

r(x) 


la  of  rlaaa  Lip  [A (6))  whenever  f  c  Lip  (A(B)l. 
Let  ua  eat  °  ° 

„  .  1  1  .  _  C-lo  .  .  t-lo 

0  ‘  SI5F  *  7S3  •  *  *  *  TIS T  • 

Under  thla  transformation  (3.21)  bacomes 

a^OljlOfjd)  b,(C)  rl  *i(7)*j(T)dl 


bj(L)  fl 

tl  *7 


T)(l-0 


-1  «  t  <  I.  where  a, (O, »>,((>,  ZjtO.TjtC)  are  related 
to  a(x),b(x),Z(x),r(x)  raapectlvely  by  a  relation  of 
the  fora 

V°  •  -l£&)  •  Z,.2t) 

f)(o  *  '(&)•«  •  *i(°  * 

Note  that  a(  and  b(  are  analytic  continuation!  Into 
»  of  i  and  b  respectively  ao  that  they  arc  In  tha 
class  LlPgl-1'1)  and  satisfy  |a?  ♦  bj|’  t  J  >  0  o« 

1—1. 1 1 .  Also,  place  f  «  lip  (2)  It  follows  from 
(3.74)  that  f(  •  Llpal-l,l).  Consequently,  In  view 
of  (3.13)  and  (3.1b),  the  solution  g.  to  (3.23)  lo  of 
claaa  Llpal-l,l). 

Thle  groves  that  Rf  la  analytic  In  Jb  and  of 


dona  Llpu(£  ) . 


Finally,  let  ua  act 


*v»> 

■a* ll ( * )  I 
(*«*) 


!(««)(*) I 


Clearly  ur  oust  have  llRU*'  *i  tor  If  •>  then 

as  a  consequence  of  the  fact  that  tha  auprsmta  on  the 
ri|ht  hand  aide  of  (3.23)  la  taken  over  a  compact 
familv,  there  would  .exist  a  ,  <  *a(J6 )  with  |g(x)|  x  1 
in  fS  such  that  Rg  4  I  (Jb),  contradicting  Kg  t  » L(S) 
for  all  g  c  >0(S).  I 

Remark  3.4.  Ve  have  also  shown  in  the  above  proof 
that  if  77  1*1-1,11  then  since  TP  •  -RQF,  we  have 
sup  |(TF)(r)|  •  oupJ(RQF)d)|  *  |>^|q8|fL.  <>•**> 

so  that 

*  l«l#M<*.  o.iT) 

where  |r|L  la  defined  In  (3.7)  —  M  la  hounded  as 

in  (3.11) . 


Map.  if  K  la  the  aaalleit  conatant  auch  that 

IfTFXa)  -  (TF)(y)  |  i  x|*-y|“  (!•»> 

(or  all  F  tilth  |)r|L  <  1,  It  folloua  that  we  bum  have 
0  <  K  <  •  and  that  for  arbitrary  F  <  L'|-l.l|. 

|  (TF)  <m>  -  (TFKy)|  «  x|Fl.U-yl"-  0.2*) 

1.2  Bpproalaot ton  on  l-l.ll. 

With 55  and  p  defined  aa  In  (2.14)  let  ua  alao 
write,  for  a  given  poaltlve  Integer  h, 

^  -  (ekh-l)/(ekhel).  k  •  -H(1)N,  1 

C-N.1  *  Ski  *  1  ;  ( 

*.(*)  *  a.  <«)/e'(r  ),  k  •  -N(I)N.  V 

k  k  k  „  7  O.JO) 

V_N_, C*>  •  (l-«>/2  -  ^(l-tk)*h(*)/2,  / 

VlU)  *  (1+,)/2  -  I  o*V*k<,)/2,  J 

Clven  F  <  l“l-l,l]  ve  alao  act 

N+l 

f„F  -  l  F<;k)bk.  (J.J1) 


Kayirh  3.V,  We  have,  with  LF  defined  by  (2.52) 

N 

F  -f  F  -  F  -  LF  -  l  0.32) 

N  k— n” 

ao  that  if  F  <  Bo(JB)  then 

IlF  -  e  Fl  -  it  h  *  >/"S>  (J  jj. 

'  "  *"  [0(NV«p(-(«<l3N),i)),  (3-33> 

1/  h  *  where  y  and  S  *  d(y)  are  poaltive 

conatant*. 

1.3  Approx lo.it ion  of  the  Integral  Equation. 

Ur  ah.ill  obtain  an  approsiaate  aplutlon  of  the 
integral  equation 


*  -  T?.  %  -  P„TF  -  V  0  3 

Then  by  Lraaa  3.1,  *  i  1  (j& )  and,  by  iaaark  1.4, 
|T|^<  ■>  ao  that  e 

au^  |v(»)|  *  |T^|rL  <  -■  (3.3 

Hence,  Belling 

u  •  v  -  Lv,  (1.4 

where  L  fa  defined  In  (2.12),  we  have 

aup  |u(t)|  i  aup  |v(t) |  ♦  aup  |(Lv)(a)|  t  4TUrL- 

*it>  *<a 


Hence,  by  laaa  2.13, 

|u(*)|  *  C|l-*2|“.  *«b.  (3.42: 

where  C  la  a  poaltive  conatant.  Indeed,  In  view  of 
(1.21),  we  nay  take  C  •  k|fL-  Combining  with  (2.21) 
we  have  for  all  a  c  (-1,1) 

|u(a)  -  I  u(*k)ak(«)/d'(«k)|  i  K|rLj/(2»alnh(»d/h)) 
^  (3.411 

Vh,r*  3  -  -up  [  ItdCifcl  <  -. 

-X«<l  ^  I* 

Furthcrnorc,  by  Uaa  2.11  and  (3.42)  with  C  *  kJf^, 
we  have 

-  2U2olt|FLv,,? 

k-IKl  (leekh)^° 

*  21+2“g|FL  ehi 

k»H+l 

,U2a  h 

«— SR2-*  «I'L-  0.45) 


f  (af )/(rZ)  -  h»(f/(rZ>)  ♦  bl^.,  =  g  (1.14)  Cooblntng  (1.45)  and  (3.41)  we  get,  for  all  a  c  (-1,1) 


ehrre  P,  .  it  i  polynow la 1  of  degree  r-1.  Let  ut 
repla<  e  (3.14)  by  the  equation 


fil  *  TH^i  "  *N 
hi  •  I  ,Vk- 


1  ckl‘k' 

k-*N-'  •  (1. 
(?N*  and  Tk  •  S^T. 


Lems  1-4.  Let  T,  T  be  conaldrrrd  as  operators  on 

rvinr  n>en 

|:T-3n|  -  aup  |(Tr)(a)  -  (T  F)(*)| 

"  -I*  1  " 

VL  ■ ' 

(o  (rapf-fh**)).  If  h  •  y/H*. 

(0  (h'Vapf-fvdaN)^) ,  (1.11) 

If  h  •  l*d/(jh)  | **,  whvrr  y  and  6  are  poaltive  niwbrr* 
auch  that  given  y  >  0,  there  ealata  6  •  4(y)  >  0. 


|u(a>  -  l  u(a.  )pk(a) | 
k*-h  * 

f  j  2l+2af-(aN-l)h J 

*  j 2*alnh(ad/h)  *  oh  j 

(1.46! 

Tliat  fa.  fron  (1.32),  (3.40)  and  (3.18), 

|v-  fkvt  -  |TF-  f,TFL 

(  J  jH2o  -(t«H-l)h) 

1  j  2*elnh(«d/h)  *  ah  j  *)rL- 

(1.47: 

Equation  (3.37)  nou  (ollowa  fro*  (3.47)  occordli^  to 
the  trier tlon  of  h.  f 

1.4.  Convorgonca  of  aeoroalnatlona  when  (1-T)~3  oaloti 

The  following  raault  lo  duo  to  Banach  (too,  for 
eaoaplo,  |7J). 

Lew*  1,7.  Lot  X  be  a  Banach  apace,  with  T  aod  (I-T)”' 
contlnuoua  lleoar  operatora  on  X.  Lot  T„  bo  a  rnwnu 


0.48) 


linear  operator  on  X  such  that 

|l  <I-TJ_I  |*ffT-THK  <  1.  0.48) 

Than  (l-T^)*1  aniata  and 

— :  •  <3‘«> 
i-Io-t)  'IIt-tJ 

Umo  3.8.  Ut  X.T  and  T_  ha  drflnrd  aa  In  Lcmi  3.7. 
Than - 

(1-T)”*  -  <t-TN)'‘ 

•  (l-T**,(T-TN)|l*0-T>',«T-TK'r,(l-T)',.O.M» 
io  that  for  *  plvrn  %  •  X,  If  (3.4fl)  t*  sil  Ixf  l«*d  thru 

il(l*T)''l  -  <1-V*'8l 

*  - !! - , -  .  <3. SI) 

Prool .  Sa»  lkcbe  (6) . 

He  can  now  prove  the  following  theorem. 

Theoren  3.4.  If  g  <  and  If  <  1-T)  mini  a  aa  a 

bounded  operator  on  L*  1-1,1)  then 

|(l-T)*'g  -  (1-1NS-1 

Ol^enpl-lndjH)'*).  h  -  (sd/.iN)'\ 

•  u  c  0.52) 

O^ni't-ilO),  h  •  y/fi4, 

where  y  and  ‘  are  positive  numbers. 

Proof_.  Sln-»  p  >  b^tjb)  we  have  g  i  1.  (-1,1 1; 
furtheraerr  J1  g  r  L"|-I,l|.  Since 

O-T)'1*  -  (I-TH)''fNg 

*  l(l-T)"‘g  -  <I-TH)-*g}  ♦  <I-TN)_,«g-  •*„*). 

O.S2)  follows  froa  0.51).  (3.49)  and  (3.33).  # 

3.5,  Calcrhln-Nyatrom  Quadrature  Approx I— l Ion . 

In  addition  to  (a),  (b)  and  (c)  of  Aaauaptlon  3.1 
lot  ua  al ao  aasuar  that 

(d)  for  rath  fixed  x  •  (-1,1),  k(x,*)  la  In  B(JD) 

(aer  definition  2.2).  and  •orcovcr  for  -1  <  t  <  1 
and  x  »  | -1.1) 

|k(a, i) |  s  Cjll-t2!0*'  (3. S3) 


where  Cj  la  a  positive  constant  and  where  a  lx  the 
ea«r  at  In  Aaauaptlon  3.1. 

Lcum.i  3_.  10.  If  F  ,  g(£> )  and  If  |F(t)|  •  C-d-t2)0”* 
on  C -1 . I )  then 

l  f  1  H4l  _  jh  I  1h  , 

f  7(t)di  -h  l  l£— —  r  r~' 

'l-i  J— N-l  (iee,h)*  lc)h*i  1 

f  0<'»P(-<«di»)S).  If  h  -  (ad/ONl)1’. 

•)  1.  L  0*> 

l  0<««P <-*«  >).  If  h  -  y/fd. 

Where  f  III  poaltlve  eonatanl  and  6  -  6(y)  >  0. 

Proof .  See  SCenger  |9,  Example  4.8),  f 

llallarly  we  can  une  (2.2S)  10  approximate  (Df)(it) 
on  (—1.1)  chla  appro.laat fon  taking  (he  lorn 


I  f*P(t)dt  -  hde*0*)  T  eJhd-(-l)J‘l>  -felh-0 

•1,^7  «  W’ 


where  the  ,ub  My  he  ralnitid  to  within  a  relative 
error  ef  the  ease  order  of  nagnltude  aa  that  oa  the 
right  hand  aida  of  (3.S4)  by  uaing  Algorithm  2.7. 

Me  shall  next  show  that  under  our  assumptions  on 
k(x.t).  our  Galerkln  approximation  schame  (3.30)  and 
(3.35)  for  aolvlng  (3.34)  aiy  ba  readily  reduced  to  the 
Myatroa  trheae  (see.  for  exsapla,  |1)  and  (6)). 

To  illuairate.  wc  observe  that  while  Interpolates 
f  at  rk,  k  •  -M(1)N,  at  well  as  at  il.  the  absolute 
value  of  the  error  Introduced  on  replacing  (u^F) (*-.,) 
bv  (f^FHI)  la 

KP.FHx^i)  -  (P„r>(l)|  -  UIK,^.1l<i-t1H1)/2 

<  lc.+l“c-M-l l«apl-(h*l)h] ,  (3.S4) 

this  being  of  aaaller  order  than  the  error  of 
Interpolation  given  by  (3.33),  which  we  have  chosen 
to  Ignore.  Similarly  *  (Pj,F)(-l)| 

la  also  bounded  by  the  right  hand  side  of  (3.S6). 

Sfallarly,  the  quadrature  formula  In  (3.S4) 
yields 

N+l 

rtf;  PyFH*)  5  h  I  2eJhk(x.t,)c./(lteJV  -  « 

T  "  j— M-l  3  3 


r  .  2l.c'3(N+,)hlk(x. *_„_,)  -  kfx.t^,))  » 

*  ,e-N.l*cMel,/U4*’t"4l>h,J-  ,3  57) 

ao  that  •  to  of  such  Mailer  order  than  the  order  of 
quadrature  error  In  (3.S4).  Hence  ue  can  Mfely  net 
<  •  0  to  get  a  much  simpler  approxlMtlon,  Which  la 
the  My  at  row  approxlMtlon. 

Similarly  if  f  <  >.(£>  then  both  u  •  af/(rZ)  and 
v  -  bX(f/(tZ)l  arc  In  k  (ft)),  nnd  the  ewalMtlon  of 
3t(f/(>2))(>)  via  Algorithm  2.7  la  straightforward  for 
x  •  but  My  be  Impossible  for  x  *  11.  Hence 

It  la  aTsoconvenlent  here  to  replace  m  •  11  by 
x  »  respectively.  The  modulus  of  the  error 

In  Mkfng  these  replacements  la  given  by 

Ivdl-vd^,)!  *  R 1 1  ~*|w.|  1°  *  2°K exp  (-o(Mfl)b), 

ao  that  it  la  of  aMller  order  than  the  Interpolation 
error  in  (3.33). 

The  exact  linear  system  corresponding  to  (3.3S)  la 
(l-T)c-t  (3.58) 

where 

£  *  tc-M-l . cIH1>T*  *  ’  . *Mel)  * 

1  here  denotes  the  unit  Mtrla  of  order  (2M+3)  and  T 
la  a  square  Mtrlx  (t  ),  n.n  .  (-M-l)(l)(H,l)  of 
order  (2N+3).  He  m5* 

S  ■  TiQfUJ  *  >«.><*  4  f«-lM  <3”> 

for  m  -  (-M-I) ( I ) (H*l) ,  and 


fl  iUJII.  .1) 

•  -  L  <7TTmO  -M-.vWlk<-.,)i(V)V,)dt 

1  ■  "  (3.60) 

(or  a.n  -  <-N-l)(l)(N+l).  and  C  i»  defined  by  (3.30). 
In  elan  of  the  above  discuaalon.  we  replace  the 
ayaten  (3.38)  -  (3.60)  by  a  perturbed  linear  ayaten 
of  the  tape  order,  given  by 


(I  -  T)  c  -  * 


where  I  la  aa  before,  e  •  (e , . . .  ,e ^l  and  the 

entrlea  g^.a  •  (-N-l)(l)(M*l)  of  tlie  vector  g  and 

t  of  the  Matrix  T,  where  n.n  •  ( — N- 1 ) ( I )(N+1) ,  are 
•,n 

given  by 

a(r  )f(r  )  . 

5.  •  TT^ziTJ  *  b<V**  '7z><V  *  •  (3-“> 

nh  (i(z)k(z.z)  *) 

*  ("iV*)'5  l  r(Vz<V  "  *  b<In*|k<*,,t.)l(v]  * 


The  singular  integrals^!  ]  in  (3.62)  and  0.63)  are 
appro* lasted  via  the  substitution  of  the  formula  on 
the  right  hand  aide  of  (3.33)  into  Algorithm  2.7. 

Bounds  on  the  difference  || c -c  ||  of  the  solutions 
of  (3.38)  and  (3.61)  can  be  obtained  by  the  well 
known  results  of  linear  algebra  see,  for  cxaaplr, 
Wilkinson  [10J,  where  we  note  that  in  tills  case  when 
h  •  (*d/(.tN))*  the  perturbations  on  the  coefficients 
of  the  aatrli  and  the  right  hand  aide  are 
0(N Vap(-(wdiN)^)). 

6 .  An  Exaaple 

Tl»e  alaorltha  described  by  (3.61)  -  (3.63)  has 
been  applied  to  the  equation 


x) 


I*  )dt  ♦  I  f*  IJlili1. 

Lj  *“*  "  JL|  t4»a-J 

•  -1  4  y\2  4  6t  4  8)“  *. 


For  this  equa:l«*n  we  h-ive  r(x)  •  I,  the  fundaaenta) 
function  2(x)  •  2(1 -a*)  »  and  the  Index  *  ■  2. 

If  w<  rhoosi-  the  particular  solution  so  that  F  ■  0.3 
at  any  two  points  of  (-1,1)  then  we  find  that 

F(x)  •  0.5  for  all  a  €  (-1,1).  By  choosing 

h  ■  x/(2N)**  and  various  values  of  W  we  have  obtained 
1  Significant  figure  In  the  appro* 1 wale  solution 
wh»»n  K  •  4,  two  significant  figures  when  W  •  8,  three 

when  N  -  16  and  five  when  h  •  32.  In  fart  we  find 

that 

■**  |F(*.  )-FM(r.  )|  •  J.0243eap{-2.2556M1,J  ,  (4.2) 

k*-N( I ) N  «na 

approxlaatrly. 
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